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Nature of the Vacuum inside the Color Flux Tube 

F. Gliozzi a * and S. Vinti a 

a Dipartimento di Fisica Teorica, Universita di Torino, 
via P. Giuria 1, 10125 Torino, Italy 

The interior of the color flux tube joining a quark pair can be probed by evaluating the correlator of pair of 
Polyakov loops in a vacuum modified by another Polyakov pair, in order to check the dual superconductivity 
conjecture, which predicts a deconfined, hot core. We also point out that at the critical point of any 3D gauge 
theories with a continuous deconfining transition the Svetitsky-Yaffe conjecture provides us with an analytic 
expression of the Polyakov correlator as a function of the location of the probe inside the flux tube. Both these 
predictions are compared with numerical results in 3D Z2 gauge model, finding complete agreement. 



1. INTRODUCTION 

The underlying assumptions of the dual super- 
conductivity (Tj of gauge theories, and its appro- 
priatenss for describing quark confinement, are 
not rigorously founded, and it is necessary to per- 
form precise numerical or analytic tests of this 
conjecture whenever possible. 

The internal structure of the color flux tube 
joining a quark pair provides an important test of 
these ideas, because it should show, as the dual 
of an Abrikosov vortex, a very peculiar property: 
it is expected to have a core of normal, hot vac- 
uum as contrasted with the surrounding medium, 
which is in the dual superconducting phase. The 
location of the core would be given by the van- 
ishing of the disorder parameter ($a/(^)) = 0, 
where $m is some effective magnetic Higgs field. 

In a pure gauge theory, the formulation of 
this property from the first principles poses some 
problems, because no local, gauge invariant, dis- 
order field &m(x) is known. As a consequence, 
one cannot define in a meaningful, precise way 
the notion of core of the dual vortex. 

A possible way out is suggested by the fact that 
in a medium in which (&m) = the quarks should 
be deconfined, then it is expected that the in- 
terquark potential inside the flux tube gets mod- 
ified. As a consequence, one may try to define a 
gauge-invariant notion of normal core of the flux 
tube as the region where the interquark interac- 



tion mimics a deconfined behavior. Of course one 
cannot speak of a true deconfinement, as it would 
require pulling infinitely apart the quarks, while 
the alleged core has a finite size. 

A simple, practical way to study in a lattice 
gauge theory the influence of the flux tube on 
the quark interaction is based on the study of the 
system of four coplanar Polyakov loops Pi, P2, P3 
and P4 following two steps 

Modify the ordinary vacuum by inserting in 
the action the pair P3 , P| acting as sources 
at a fixed distance R. 

Evaluate in this modified vacuum the correla- 
tor the other pair Pi , P| of Polyakov loops 
which are used as probes. 

The correlators in the two vacua are related by 



(PiPjPaPj) 



(1) 
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In this note we study some general properties of 
these correlators for T < T c . In particular, we 
point out that at T — T c the functional form 
of these correlators is universal and in some 3D 
gauge theories can be written explicitly, even in 
finite volumes. 

2. FOUR POLYAKOV LOOPS 

Consider the system of four parallel, coplanar 
Polyakov loops, symmetrically disposed with re- 
spect the origin of a cubic lattice with periodic 
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boundary conditions in the direction of the imag- 
inary time (which coincides with the common di- 
rection of the loops). We study their correlator 



(PiPtPsP]) = (P(- 



f)^(|)> (2) 



as a function of r < P. For large R and r ~ R it 
obeys the asymptotic factorization condition 



(3) 



When T < T c , assuming the usual area law 
{P\P \) oc exp(— ar/T), where a is the string ten- 



sion, yields 

(PPj)^ ~ exp(ar/T) ~ l/(P 1 P 2 t > 



(4) 



which gives an apparent repulsion between the 
two probes due to the attraction of the two 
sources. 

The other limit r <C R is more interesting, be- 
cause the kinematics does not force any factoriza- 
tion and different confinement models suggest dif- 
ferent behaviors. In particular in the naive string 
picture one is tempted to assume the factorization 
(^) even in this limit, because within this assump- 
tion the total area of the surfaces connecting the 
Polyakov loops is minimal. On the contrary, in 
the dual superconductivity it is expected that the 
test particles probe the short distance properties 
of the hot core of the flux tube, thus the correla- 
tor in the modified vacuum would approach to a 
constant 



(P)t>t ) from above and 



<PiP 2 %- > (PiPl) 



(r<P,T<T c 



(5) 



In the range T > T c the interior of the flux tube 
is in the same phase of the surrounding region 
and the mutual interaction between the two near 
probes should not depend on the presence of very 
far sources, then 



(PiP 
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(PiPj) (r«P,T>T c 



(6) 



2.1. Critical Behavior 

According to the widely tested Svetitsky-Yaffe 
conjecture, any gauge theory in d + 1 dimensions 
with a continuous deconfining transition belongs 
to the same universality class of a d-dimensional 
C(G)-symmetric spin model, where C(G) is the 
center of the gauge group. It follows that at 



the critical point all the critical indices describing 
the two transitions and all the adimensional ra- 
tios of correlation functions of corresponding ob- 
servables in the two theories should coincide. In 
particular, since the order parameter the gauge 
theory is obviously mapped in the corresponding 
one of the spin model, the correlation functions 
among Polyakov loops should be proportional to 
the corresponding correlators of spin operators: 

(Pi ... P2n)T=T c OC (si ... S 2n ) ■ (7) 

Conformal field theory has been very successful 
in determining the exact form of these universal 
functions for d = 2 even in a finite box, which is 
a precious information for a correct comparison 
with numerical simulations. 

In particular, using the known results of the 2D 
critical Ising model in a rectangle L\ x L 2 with 
periodic boundary conditions we can write ex- 
plicitly the correlator of any (even) number In of 
Polyakov loops of any 2+1 gauge theory with 
C{G) = 7L 2 . Let Xj,yj be the spatial coor- 
dinates of Pj and define the complex variables 
iL 2 /L 1 . Then 



z,; = r 2 - + iy 2 - and r 



(Pi...P 2 „) 2 = C „^ J2 ( 8 ) 



E 

U=l Ej=±l 



with e ■ z = e i z i an d the primed sum is con- 
strained by J2i £j = ; c„ is an overall constant 
that can be expressed by factorization in terms 
of c\ . The universal functions A v and P^ can be 
written in terms of the four Jacobi theta functions 
i?„(z, t) as follows 

aej/2 

(9) 



B 



■&l{Zi - Zj,T) 



i?i(o,T) 



0„(O,t) 

In the infinite box limit L%, L 2 
Taylor expansion 

$ v (z, t) = a v {\ - 6 ljV ) + b„z + 0(z 2 ) 
the correlator (0) becomes 



(10) 

co, using the 



(PiP 2 P 3 P 4 ) = 



4c? 



R 



(Rr)i V R-r 
which satisfies both factorizations (pph. 



(11) 



(12) 
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3. CLUSTER ALGORITHM 

In order to test the above formulae at critical- 
ity it is convenient to perform the numerical sim- 
ulations in the simplest model belonging to the 
above-mentioned universality class, which is the 
the 3D 7L2 gauge model. 

Using the duality transformation it is possible 
to build up a one-to-one mapping of physical ob- 
servables of the gauge system into the correspond- 
ing spin quantities. A great advantage of this 
method is that it can be used a non local cluster 
updating algorithm ||, which has been proven 
very successful in fighting critical slowing down. 

In this framework it is easily shown that the 
vacuum expectation value of any set {C\ . . . C' n } 
of Polyakov or Wilson loops of arbitrary shapes 
is simply encoded in the topology of Fortuin- 
Kasteleyn (FK) clusters: to each Montecarlo con- 
figuration we assign a weight 1 whenever there 
is no FK cluster topologically linked to any 
Ci G {Ci . . .C n }, otherwise we assign a weight 
0. Let No and Ni be the number configurations 
of weight and 1 respectively, then we have sim- 
ply 

Ni 



(Ci ■ • • C„> 



No + Ni 



(13) 



This method provides us with a handy, very pow- 
erful tool to estimate the correlator of any set of 
Wilson or Polyakov loops even at criticality. 

4. RESULTS 

In order to test the critical behavior of the 
multiloop correlator one has to know with high 
precision the location of the critical tempera- 
ture as a function of the coupling (3. We took 
advantage of ref. where these critical values 
have been obtained with an extremely high sta- 
tistical accuracy. We report in Fig.l some re- 
sults at (3 = 0.746035 corresponding to l/aT c = 
N tc = 6 and to a string tension aa 2 — 0.0189(2). 
The open circles are the data for the correlator 
(P(-§)P(§)) in a N t x N x x N y lattice with N t = 
3N tc ,N x = N y = 64. They are well fitted by 
the one-parameter formula cexp(— aN t r)/rj(i^-), 
where the Dedekind 77 function takes into account 
the quantum contribution of the flux tube vi- 



brations §. The square symbols correspond to 
the correlator (P(— §)-P(|))gg at the same tem- 
perature, in presence of a pair of sources at a 
distance R = 24a. The data in the central re- 
gion are well fitted by a two-parameter formula 



c qq exp(-aN t r)/r)(iq±) + b 



b which 



•^qq^fy ji vv 2r > 1 " ' ^qq 

simulates a high temperature behavior with a 
screening mass m = aN t and an order param- 
eter (P) = \fb. The black circles correspond to 
(P(-|)P(|)P(-f)P(f )) evaluated at T = T c 
with R — 16a. They fit nicely eq.(||) (continu- 
ous line). Note that such a curve does not con- 
tain any free parameters, being C2 = V^2c| with 
1 

c\Nx — 0.199(4) as estimated by measuring 
(P1P2) on lattices of different sizes at T = T c 
and N tc — 6. 
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Figure. 1. Correlator of two Polyakov loops 
inside and outside the flux tube 
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